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. Abstract. Magnetically mediated disk outflows are a leading paradigm to explain 

' winds and jets in a variety of astrophysical sources, but where do the fields come 

£SJ \ from? Since accretion of mean magnetic flux may be disfavored in a thin turbulent 

■ disk, and only fields generated with sufficiently large scale can escape before being 

shredded by turbulence, in situ field production is desirable. Nonlinear helical inverse 
dynamo theory can provide the desired fields for coronae and outflows. We discuss 
the implications for contemporary protostellar disks, where the MRI (magneto- 
rotational instability) can drive turbulence in the inner regions, and primordial 
^s^j , protostellar disks, where gravitational instability drives the turbulence. We empha- 

^ • size that helical dynamos are compatible with the magneto-rotational instability, 

and clarify the relationship between the two. 

(N ' 
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1. Introduction: the dynamo-corona-outflow paradigm 



5_i ■ Bipolar jet-like outflows are commonly observed in protostellar sys- 

"to . terns (Richer et al. 2000), as is enhanced X-ray activity (Feigelson & 

Montmerle 1999). Though the correlation between outflows and X-ray 
activity is not entirely systematic, magnetic fields probably play an 
^ | important role for both phenomena. Large scale outflows are likely 

5_i ' magneto-centrifugally "fling" driven (Blandford &: Payne 1982) or mag- 

netically "spring" driven (e.g. Lynden-Bell 1996) from large scale open 
field lines, while coronal dissipation is likely the result of reconnection 
from closed loops, perhaps as they open (relax) to form larger scale 
structures. Both jet and coronal fields are plausibly the result of large 
scale fields produced inside protostellar disks because (1) magnetic flux 
is difficult to accrete in a thin turbulent disk (2) only large scale fields 
are able to escape from the disk before being shredded by turbulence 
therein. A working paradigm is this: helical dynamo — > large scale 
fields — > coronal fields — > field lines that relax and open up — > magne- 
tocentrifugally launched outflow with a significant supersonic vertical 
velocity component — > asymptotic outflow opening angle whose tangent 
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is the ratio of the expansion speed to the vertical wind speed. The angle 
may be further reduced by magnetic collimation. 

We focus here on understanding the principles that govern the strength 
of large scale disk fields produced by dynamos, and on the mechanical 
luminosities of the resulting winds. We do not focus on the possible 
role of stellar dynamos or the accretion ejection instability (Tagger & 
Pellat 1999; Varniere & Tagger 2002). 

2. Need for in situ generation of large scale fields 

2.1. Small scale fields shred before they escape 

For magnetic fields to launch jets or power coronae, the field buoyancy 
time tb, must be < its diffusion time in the disk. Since tj, <i h/Uj,, 
where h is the disk 1/2 thickness and Ub is the buoyancy speed (< 
Alfven speed U a) associated with a structure whose smallest gradient 
scale is L, we require h/Uf, < L 2 /(3 for escape, where (3 ~ U2I is the 
turbulent diffusion coefficient for the magnetic field, and U2 and I are a 
measure of the dominant turbulent speed and scale, where (e.g. I ~cube 
root of the turbulent cell volume {l x lyh) 1 ^)- The escape condition can 
be rewritten L 2 > Ih(u2/Ub) and since the ratio on the right (more on 
this later) is <: 112 /U a ~ 1, the condition L 2 > Ih is required for escape. 

Now consider a Shakura-Sunyaev (Shakura & Sunyaev 1973) disk 
viscosity v ss ~ (3 = U2I ~ a ss c s h ~ u?,/Cl ~ u A2 /^, where ua,2 is 
the Alfven speed associated with the turbulent -B-field and U2 ~ ua,2 
follows from simulations (Hawley et al. 1995; Brandenburg et al. 1995). 

Then a ss ~ u 2 A2 /c 2 s or ua,2 ~ all 2 c s , so that I = aj^h. The escape 

condition L 2 > Ih then gives L > all 4 h. Note that the minimum L 
applies to the smallest dimension of a given structure since this deter- 
mines the diffusion time. This escape condition is not easily satisfied 
in the nonhelical MRI dynamo (see minimum k y ,k z in fig. 3a (from 
Hawley et al. 1995)). More on this point later. 

2.2. Mean flux is not frozen: diffusion may beat accretion 

Can the large scale fields required by jets simply be accreted? Though 
3-D MHD simulations are needed, it may be difficult for thin disks. 
The reason is that turbulent magnetic diffusivity [3 ~ may keep 
the large scale field from being frozen on an advection time scale. (In 
general, the extent of flux freezing even for the total field depends on 
the turbulent spectrum. Turbulence produces many small scale stru- 
cures which thus enhances the importance of the diffusion term in 
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Figure 1. Saturated dynamo energy spectra (Maron & Blackman 2002). 

the induction equation.) To see the potential problem, we compare 
the ratio of the rate of flux advection to that of diffusion. Call d the 
dominant vertical variation scale of the field strength. The dominant 
variation scale of the radial velocity is r. The above mentioned ratio 
is then V x ( v x B)/(/3V 2 B) ~ (d 2 /r 2 )(a ss hc s / (3). The diffusion term 
thus dominates when the plausble (though not proven) relations d < r 
and a ss hc s ~ (3 apply. That diffusion can be faster than advection is 
supported by numerical 2-D calculation in (Lubow et al. 1994), but 
future 3-D numerical testing is required. This motivates an in situ 
helical dynamo, whose growth rate can exceed the diffusion rate. 



3. MHD In Situ Amplification: Direct vs. Inverse Dynamos 

We classify dynamos into two types: direct and inverse. Direct dy- 
namo action sustains magnetic fields on scales at or below the dominant 
turbulent scale. This does not require helical turbulence although helic- 
ity does influence the magnetic spectrum. In contrast, inverse dynamo 
action describes amplification on scales larger than that of the dominant 
scale of the turbulence. The label "inverse" is used to suggest an inverse 
cascade. The turbulence must be helical to generate and sustain large 
scale flux over times longer than an eddy turnover time which can in 
turn escape to coronae and drive outflows. The distinction is illustrated 
in Fig 1. (Maron & Blackman 2002) for forced turbulence. The mag- 
netic and kinetic energy spectra are plotted for different values of the 
fractional helicity fh = |u2 • Vxu2 \/\u?,k2 1. The nonhelical forcing at 
&2 = 4.5 produces no magnetic energy at wavenumbers k < &2, whereas 
the fh = 1 case produces a large peak at k = 1. 
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3.1. Direct dynamo 

For nonhelical direct dynamo action, the field grows by random walk, 
field line stretching (Kazantsev 1968; Parker 1979). The turbulence can 
be externally driven by isotropic forcing (e.g. Maron & Cowley 2002; 
Haugen et al. 2003) , or self-generated by an angular velocity gradi- 
ent (e.g. Balbus & Hawley 1998). In either case, turbulent stretching 
compensates for exponential decay from turbulent diffusion. The latter 
operates because the random motions of the gas also mix the field 
lines, inducing a cascade to small scales where dissipation occurs. In 
3-D, a steady state balance can be achieved. In the saturated state, 
the magnetic energy ~ turbulent kinetic energy integrated below the 
dominant turbulent scale. How far the magnetic energy peak is from the 
forcing scale and the overall nonhelical dynamo spectrum is an active 
area of research (e.g. Maron & Cowley 2002; Schekochihin et al. 2002; 
Haugen et al. 2003), which we do not discuss further here. 

3.2. Helical inverse dynamo 

The helical inverse dynamo amplifies field on scales larger than that 
of the turbulence. This is most desirable for coronae and outflows. Fig 
2a. shows a traditional a£l kinematic dynamo diagram (Parker 1979). 
Consider an initially weak toroidal (=encircling the rotation axis) loop 
of the magnetic field embedded in the rotator. With a vertically decreas- 
ing density gradient, a rising turbulent eddy threaded by a magnetic 
field will twist oppositely to the underlying global rotation to conserve 
angular momentum. Statistically, northern (southern) eddies twist field 
clockwise (counterclockwise). This is the "a^" effect and the result is 
a large scale poloidal field loop. Differential rotation shears this loop 
(the "fT-effect). The bottom part reinforces the initial toroidal field 
and the top part diffuses away. The result is exponential growth. The 
ad effect can be supplied by any source of 3-D turbulence in a stratified 
rotator, including gravitational instability or the MRI. 

This dynamo is revealed by averaging the magnetic induction equa- 
tion over a local volume and breaking all quantities (velocity U, mag- 
netic field B in Alfven velocity units, and normalized current density 
J = VxB) into their mean (indicated by an overbar) and fluctuating 
(lower case) components. The result is 

9,B = Vx(a d B + UxB) + (^ + A)V 2 B, (1) 

where A is the magnetic diffusivity. The U term incorporates the f2- 
effect, the (3 term incorporates the turbulent diffusion (we assume 
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Figure 2. Kinematic helical dynamo diagram (left) ignores magnetic helicity con- 
servation. In dynamic theory field should be replaced by ribbon, showing the 
conservation (Blackman & Brandenburg 2003). 



constant (3) and the first term on the right incorporates the a^-effect. 
Accretion disk a— Q dynamo models have been proposed (Pudritz 1981; 
Reyes- Ruiz & Stepinksi 1997; 1999; Riidiger et al. 1995; Brandenburg 
& Donner 1997; Campbell 1999; Riidiger & Pipin 2000) but without a 
fully dynamic theory of the field saturation value. We now discuss this 
for forced turbulence, and expand further for disks in section 4.2. 

For forced helical turbulence, a<2 is (Blackman & Field 2002) = 
CKdo + c^m where a^o = — (r/3)u • Vxu, and a m = (r/3)b • Vxb, where 
r is a turbulent damping time. The a m correction term arises from in- 
cluding the magnetic forces that backreact on the velocity (Pouquet et 
al. 1976; Blackman & Field 2002). From the form of a^, it is evident that 
growth of a m can offset a^o and quench the dynamo (Field h Black- 
man 2002; Blackman & Brandenburg 2002; Blackman & Field 2002). 
Magnetic helicity conservation determines the dynamical suppression. 
The magnetic helicity, a volume integral H = J A ■ TidV = (A • B)V 
satisfies (Woltjer 1958; Berger & Field 1984) 

dfH = —2\C — surface terms, (2) 

where V x A = B, and C = (J • B)V. Without non-diffusive surface 
terms H is typically well conserved in astrophysics. Since H is a mea- 
sure of field line "linkage" and "twist," its conservation implies that 
the a^-effect cannot produce a net magnetic twist, just positive and 
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negative magnetic twists on different scales. The conservation equations 
like (2) for Hi and H2 (where 1(2) indicates large (small) scales) are 

d t Hi = 2S - 2\k\Hi - surface terms (3) 

d t H 2 = -2S - 2\k\H 2 - surface terms, (4) 

where S = (a^ki/e — (3kf)Hi and a<i = cxq + a m , where e < 1 is the 
fraction of the large scale field energy which is force-free. The solution 
without surface terms and with e = 1 provides an estimate for the 
fully helical field large scale field component, which can be further 
augmented by differential rotation. Because a m oc k 2 H 2 , the solution 
(Field & Blackman 2002; Blackman & Brandenburg 2002; Blackman & 
Field 2002) shows that for initially small H 2 but large ao, Hi grows. 
Growth of Hi implies the oppositely signed growth of H 2 . This H 2 
backreacts on ao, ultimately quenching ad and the dynamo. 

Fig. 2a does not capture magnetic helicity conservation because a net 
helicity is generated. The resolution of this problem is fig. 2b (Blackman 
& Brandenburg 2003), where the field is treated as a ribbon. The large 
scale writhe corresponds to to the large scale magnetic helicity and to 
mean field dynamo action, while the small scale twist of opposite sense 
along the ribbon is the result of magnetic helicity conservation. The 
latter twist produces a force which resists the bending of the field, and 
thus encapsulates the backreaction and growth of ot m . Before discussing 
the saturation, more on the relation between the MRI and the direct 
vs. inverse helical dynamos is addressed. 



4. MRI is compatible with Direct and Inverse Dynamos 

4.1. MRI IN NONSTRATIFIED DISK: DIRECT DYNAMO ONLY 

In a sufficiently ionized disk with a weak magnetic field (note: field must 
be strong enough for growth to beat diffusion at critical wavenumbcr 
(Balbus &; Papaloizou 1999)) the MRI results from an outward de- 
creasing angular velocity gradient. Without vertical stratification, and 
thus without any pseudoscalar helicity, the MRI does not produce mag- 
netic fields on scales larger than that of the developed turbulence. The 
marginally unstable (and dominant) growth mode of the MRI satisfies 
kuA ~ ^ (Balbus & Hawley 1998) so that as the magnetic field is 
amplified from the instability, larger and larger scales become unstable 
to growth. The largest scale that can grow is that associated with the 
scale height h. Without stratification or helicity, the largest scale of the 
turbulent magnetic field matches the largest turbulent velocity scale 
(Fig. 3a (Hawley et al. 1995); k y is the toroidal wavenumber.) 
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The flux associated with the field on the largest scale in the nonstrat- 
ified MRI would change every turnover time as the field is shredded. 
The field is large scale only in the sense that the box height limits the 
scale of the turbulence NOT because there is any sustained ordered flux 
on scales larger than the turbulence. Furthermore, for a periodic box, 
the net flux cannot grow, as it is conserved. Note also from Fig. 3a, 
that the largest field structure in the vertical (k z )and radial (k x ) direc- 
tions are smaller than in the toroidal direction (k y ). This is important 
because the constraints of section 2.1 apply to the smallest dimension 
of the largest scale structures. In sum, a nonstratified nonhelical MRI 
does not produce the fields required by jets/coronae. 

Nonstratifed, nonhelical MRI driven turbulence and nonhelically 
forced turbulence behave qualitatively similarly with respect to the ab- 
sence of fields on scales larger than the dominant scale of the turbulent 
motions. Note however, that the nonhelical MRI, unlike the nonhelical 
forced case, more nearly matches the turbulent kinetic energy spectrum 
on larger scales (compare the first two panels in fig. 3 to the curves 
in fig. 1 with fh = 0). This is likely due to the fact that the global 
shear growth time is the same on all scales in a shearing box, and thus 
more able to beat diffusion on the larger scales where it can provide 
additional magnetic energy compared to the case without global shear. 

4.2. MRI IN STRATIFIED DISK: BOTH DIRECT & INVERSE DYNAMOS 

Unlike the nonstratified case, when the MRI operates in a vertically 
stratified medium it can have a built in inverse helical dynamo: the 
stratification means that turbulence produced by the MRI can be he- 
lical, in the manner discussed in section 3.2 and also from magnetic 
buoyancy (Brandenburg 1998; Riidiger & Pipin 2000). One can then 
use a nonlinear mean field theory to model the growth and saturation 
of the sustained large scale flux. This is independent of the fact that the 
MRI will always also produce direct dynamo action that amplifies fields 
as described. That is, when the MRI operates in a stratified medium, 
one gets BOTH direct and inverse helical dynamo action. 

Mean field dynamo growth is seen in magnetorotational instability 
simulations (Brandenburg et al. 1995, Brandenburg & Donner 1997). 
Fig. 3b shows a dynamo cycle in which mean flux is maintained over 
30 orbits. In that simulation, large and small scale magnetic helicities 
are seen to have opposite sign (as expected from discussion of section 
3) but (u • Vxu) has the opposite sign to (j • b) which is not expected if 
ctd = «(20+QW There is an additional buoyancy term a& ~ — | \g z (^b x )\ 
(g z is vertical gravity constant) which can exceed ado and carry the 
opposite sign (Brandenburg 1998; Riidiger & Pipin 2000). Then, for 
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example, in the northern hemisphere ay ~ —T\g z ((5p/p)b x )\ + |r(v • 
Vxv)|/3 + r(j • b)/3. The 1st and 2nd terms on the right have opposite 
signs and the sign of the large scale magnetic helicity is determined 
by whichever of those two initially dominates (since ad is the growth 
coefficient). This then predicts the sign of the last term on the right, 
which grows with opposite sign to the large scale magnetic helicity, 
eventually quenching the larger of the 1st two terms on the right. 

The fact that the stratfied MRI leads to both helical dynamo action 
and direct dynamo action, whereas the nonstratified MRI leads only to 
direct dynamo action is the key to clarifying previous discussions of the 
MRI and the tradititional mean field dynamo: It is sometimes stated 
that even a nonstratified MRI generates "large scale" fields without 
the need for helical turbulence. But this statement relates directly to 
section 4.1. The fields referred to as "large scale" in this nonhelical 
context are just the result of direct (not helical) dynamo action: 1) they 
only reach the maximum scale of the turbulence, not larger, and 2) they 
change the sign of their flux every eddy turnover time (every rotation 
period). That such fields are generated without helicity, does not make 
the MRI inconsistent with helical dynamos: as discussed above, when 
stratified turbulence is presence both direct dynamo and inverse helical 
dynamos operate together and large scale fields on scales larger than 
that of the dominant turbulence arise. 

Like Brandenburg et al. (1995), Stone et al. (1996), studied the 
stratified MRI. But the latter used vertically periodic boxes so no net 
flux could grow. Flux could have appeared in each ^ of the Stone et al. 
(1996) box with opposite sign, but this was not a large region, given 
the resolution: Measuring finite pseudoscalars in a periodic stratified 
box requires averaging only over ^ of the vertical thickness, since the 
sign changes across the midplane. For a periodic stratified box, the 
pseudoscalars must vanish at the midplane and at the top and bottom. 

5. Helical Dynamo Saturation in Disks & Outflow Power 

The helical dynamo grows helical large scale field first "kinematically" 
(independent of dissipation and with growth rate not too different from 
ado) until \a m \ ~ \a\, — ado\, after which the field growth is resistively 
limited. Setting \a m \ ~ \a^ — ado\, and using magnetic helicity conser- 
vation, then provides a reasonable estimate of the helical field by the 
end of the kinematic phase. Growth beyond that value is often too slow 
to be of interest. The large scale helical field at the end of the kinematic 
phase is (Blackman & Field 2002; Tan & Blackman 2003) 

B H ~ (4vrpA) 1 /2 M2 (/// l )V2 „ (4 7r p) 1 /2 aMCaj ( 5 ) 
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Figure 3. (a) Nonhelical MRI spectra from Hawley et al. (1995) and (b) MRI driven 
helical dynamo cycle in stratified non-periodic box (Brandenburg & Donner 1997). 



where fh is the fractional kinetic helicity, I is the dominant turbulent 

scale, h is the disk height (=the vertical scale of the growing mean 

1 11 

field), and the latter similarity follows for disks since l/h ~ a ss and 

1/2 

U2 ~ a ss c s , which in turn follow from the scalings of section 2.1. 

Eqn. (5) is an estimate of the helical field but the nonhelical toroidal 
field is further amplified by shear. The toroidal field is linearly stretched 
above B r during a buoyancy time scale. The buoyancy loss time is 
Tb ~ h{Anp) 1 / 2 / B^, so linear growth of above (5) in a time r& gives 

B^ = B r riT b = a^ 2 c s (4vr / 9) 1 / 2 so ~B r [B$ = all 2 . 

The above mean fields have their smallest variation scale ~ h and 
thus satisfy the conditions of section 2.1 and can rise to the corona, 
open up, and launch outflows. By using V • B one can use the disk 
field to estimate the field at the surface and the Poynting flux there 
(Tan & Blackman 2003). This provides an estimate for the magnetic 
luminosity available for magneto-centrifugal (Blandford & Payne 1982) 
or magnetic "spring" (Lynden-Bell 1996) driven outflows 

L mag ~^/(ExB).dS~ ^ 

~ 990 T I a ss \ 2 m* m* I 10 cm \ 
_ ZZU Lq yjQ^I J M^ 10 --2 MQ / yT { r ), 

using normalizations for primordial protostars (Tan & McKee 2003). 



6. Outflows in Contemporary & Primordial Protostars 

The time scale for which the helical dynamo operates at a given radius is 
not the viscous time, but the accretion disk's lifetime (~ star formation 
time). This is because the mean field is NOT frozen into the flow. The 
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Figure 4- (a) Dynamo saturation regime for primordial protostellar disks (b) Evolu- 
tion of radiative and wind mechanical luminosities from protostar and disk; outflow 
momentum; effect of wind on star formation efficiency (Tan & Blackman 2003). 



amplification is statistical, lasting as long as turbulence is supplied. For 
contemporary star formation, we expect that growth to the kinematic 
saturation value occurs in a small fraction of the star formation time, 
since the initial fields after collapse from molecular clouds are already 
quite strong and the star formation time scales are relatively long. Con- 
temporary protostellar disks should be able to easily sustain outflows 
powered by fields from a helical dynamo, with the helical turbulence in 
the inner regions of the disk supplied by the MRI. 

For primordial protostellar disks, although magnetic Reynolds num- 
bers are high, the seed fields are weaker and the MRI may be unfavor- 
able for the maximally unstable mode. Instead, gravitational instability 
(Gammie 2001, Wada et al. 2002; though these are 2-D simulations) 
may supply the 3-D helical turbulence for much of the dynamo evo- 
lution. In addition, we expect the star formation time scales to be 
relatively short since the accretion rates are ~ lO _2 M /yr (Tan k, 
McKee 2003). Gravitational instability in the disk may maintain the 
Toomre parameter Q $ 1 (Gammie 2001; Tan &: Blackman 2003): 
Where the disk is initially stable (such that viscosity is small), mass 
accumulates, triggering the instability. This can induce turbulent ac- 
cretion rather than global angular momentum transport (e.g. Gammie 
2001). Typical accretion rates, while large by contemporary standards, 
are small enough that fragmentation is not expected (Tan & Blackman 
2003). 

Tan &: Blackman (2003) used a primordial protostellar model (Tan 
& McKee 2003) to determine the m* vs. rii* regime for which a helical 
dynamo reaches the end of the kinematic regime (discussed above). In 
Fig. 4a., the rising curves mark the condition where the helical field 
has been able to saturate, as the disk is old enough to have allowed 
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sufficient dynamo amplification. The disk dynamo is assumed to start 
operating once the centrifugal radius of the accretion flow is ~ twice 
the stellar radius, which happens when m* init = 

0.4M© in the fiducial 

case. The dashed, solid, and dotted lines show the saturation condition 
for a ss = 1CT 3 , 10 -2 , 0.1, all with radius r = 10 13 cm and seed field of 
10~ 16 G. With a ss = 0.01, the long-dashed line shows the condition at 
r = 10 14 cm, the dashed-long-dashed line shows it for a seed field of 
10~ 26 G, and the dot-long-dashed line shows it evaluated with m*^ t = 
4.7M , corresponding to a more slowly rotating initial gas core. 

Protostellar growth starts from low masses with accretion rates that 
decline with time (and so also with protostellar mass) , as shown by the 
straight thick solid line. In the fiducial case the field in the inner disk 
saturates when m* ~ 1M & . This does not depend sensitively on a ss . In 
lower angular momentum cores the disk emerges later and the dynamo 
then rapidly saturates. The main message of Fig 4a is that dynamo 
saturation and thus dynamo-driven winds may occur at relatively low 
protostellar masses, even in the primordial case. 

Fig. 4b (Tan & Blackman 2003) shows the evolution of protostellar 
luminosities (radiative and mechanical), the cumulative specific outflow 
momentum relative to stellar mass (assuming that the wind speed is 
the escape speed from the protostellar surface and that the flow kinetic 
energy dominates far from the star), and the resulting star formation 
efficiency due to erosion of the core by the outflow for the fiducial 
case of Fig 4a, with q ss = 0.01. The efficiency calculation is based on 
Matzner & McKee (1999,2000). The initial core's axisymmetric angular 
density profile is specified by dM/dft = (l/47r)Q(/x)M, with ^ = cos# 
and Q{n) = (1 - fJ 2 ) 71 / Jq(1 — ^ 2 ) n dfi. The solid line is n = (isotropic 
core), dotted is n = 1, dashed is n = 2, long-dashed is n = 3, dot-dashed 
is n = 4. The message from fig. 4b is that dynamo driven winds can 
strongly affect the star formation efficiency for masses m* > 100M Q , 
at which the protostar has contracted to the main sequence radius, and 
the contraction has enhanced the accretion and wind luminosities. 
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